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The Design Analysis Brach (NE-M1) at the Kennedy Space Center has not had the 
ability to accurately couple Rigid Body Dynamics (RBD) and Computational Fluid 
Dynamics (CFD). OVERFLOW-D is a flow solver that has been developed by NASA to have 
the capability to analyze and simulate dynamic motions with up to six Degrees of Freedom 
(6-DOF). Two simulations were prepared over the course of the internship to demonstrate 6- 
DOF motion of rigid bodies under aerodynamic loading. The geometries in the simulations 
were based on a conceptual Space Launch System (SLS). The first simulation that was 
prepared and computed was the motion of a Solid Rocket Booster (SRB) as it separates from 
its core stage. To reduce computational time during the development of the simulation, only 
half of the physical domain with respect to the symmetry plane was simulated. Then a full 
solution was prepared and computed. The second simulation was a model of the SLS as it 
departs from a launch pad under a 20 knot crosswind. This simulation was reduced to Two 
Dimensions (2D) to reduce both preparation and computation time. By allowing 2-DOF for 
translations and 1-DOF for rotation, the simulation predicted unrealistic rotation. The 
simulation was then constrained to only allow translations. 
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Resultant force components applied at CM in x, y, and z directions respectively, lbf 
Height from sea-level, ft 

Moments of inertia about the x, y, and z axis respectively, slug • ft 2 
Products of inertia about the x, y, and z axis respectively, slug * ft 2 
mass, slug 

Y component of resultant moment applied at CM, lbf ft 
Free stream Mach number 

Engineering and Technology Directorate-Mechanical Division-Design and Analysis Branch 
Freestream dynamic pressure, psf 
Freestream Reynolds number, 1/ft (per foot) 

Second 

Freestream temperature, R 

Freestream velocity, ft/s 

Angle-o f-attack, deg 

Angle-of-sideslip, deg 

Flight path angle, deg 

Freestream dynamic viscosity, lb • s /ft 2 

Time step, s 


Computational fluid dynamics 
Chimera Grid Tools 


1 USRP Intern, Engineering Design Analysis Branch, Kennedy Space Center, Undergraduate at Purdue University, 
School of Aeronautics & Astronautics Engineering. 

2 Fluids Systems Lead, Engineering Design Analysis Brach, FL 32899, Kennedy Space Center. 


NASA Kennedy Space Center 


l 


August 12, 201 1 


NASA USRP Internship Final Report 


CM 

DOF 

GUI 

LCC 

LPVD 

NASA 

RBD 

SLS 

SRB 

2D/3D 

XML 


Center of Mass 

Degrees of freedom 

Graphical User Interface 

Launch Commit Criteria 

Launch Pad Vehicle Drift 

National Aeronautics and Space Administration 

Rigid Body Dynamics 

Space Launch System 

Solid Rocket Booster 

Two dimensional/three dimensional 

Extensible Markup Language 


A - 


314 . 5 ' 


0 

33 . 0 ' 


w— 


209 . 0 ’ 


i76.r 


1. Introduction 

T HE Engineering Design Analysis Branch (NE-M1) in the Mechanical division of the Engineering and 
Technology Directorate at the National Aeronautics and Space Administration (NASA) Kennedy Space Center 
(KSC) has never had the ability to couple rigid body dynamics (RBD) and 
computational fluid dynamics (CFD). Technical management within NE-M1 
can identify several analyses that may benefit by coupling dynamic motion and 
CFD. The requirements of this project were to demonstrate to NE-M1 the 
capability of analyzing moving bodies subjected to aerodynamic loads. 

With the growing possibility of the new Space Launch System (SLS), NE- 
M1 expects a great deal of analysis work in the future that corresponds to the 
SLS. Thus two different moving body simulations were developed based on 
geometry derived from a conceptual SLS that is being considered. Figure 1 
shows the geometry of a possible SLS configuration. The first study modeled 
was the RBD of a Solid Rocket Booster (SRB) during its ejection from the core 
stage. The second was a simulation of the SLS as it departs from a launch pad 
under the influence of a crosswind. The outcome of the two simulations in this 
project was to demonstrate to technical management in NE-M1 that CFD can be 
coupled to RBD. 

OVERFLOW 1 is the name of the software that was chosen to demonstrate 
moving body capability. OVERFLOW is a flow solver that was developed by 
NASA and is non-commercial. Thus, the code is readily available within 
NASA. Since OVERFLOW makes use of overset structured grids to compute 
fluid flow around geometries, a description and nomenclature of the overset 
gridding method will be discussed. 

The overset or Chimera grid method is a discretization technique that allows 
a physical domain to be decomposed into multiple overlapping grids. 

OVERFLOW was designed to make use of overset structured grids and thus 
any subsequent reference to grid will refer to such grids. Chimera Grid Tools 2 
(CGT) is the software that was developed by NASA Ames and was the software 
used in this project to construct all grids. 
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Figure 1. Possible SLS 
configuration. The diameters 
and lengths used to construct 
the grids for SLS concept 
vehicle. 


II. Geometry Descriptions and Modeling 

The dimensions of the SLS that was modeled for both the SRB separation and the launch pad vehicle drift 
(LPVD) simulations were derived from figure 1. Only lengths and diameters were provided and thus, specific 
contours of the payload faring and nose cones were approximated to aesthetically represent figure 1. Since the 
purpose of the simulations presented here was to demonstrate the capability of RBD coupled to CFD, the exact 
geometry was not crucial. 

As a preliminary demonstration of CFD and RBD, major simplifications were made to the model. In this project, 
rocket plumes, nozzles, and attach hardware were not included and thus any phenomena caused by these elements 
were not represented. 

Table 1 is a summary of fight conditions for the SRB separation and the LPVD simulations respectively. The 
flight conditions for the SRB separation case were derived from references 3 and 4. For the LPVD simulation, the 
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flight conditions were chosen such that the most unfavorable wind 
conditions that could occur for the SLS would be modeled at sea 
level. Based on a Launch Commit Criteria (LCC) for wind, a 
crosswind of 20 knots was simulated. 


A. Geometric Modeling of SRB Separation Simulation 

Feature curves were identified as the first step in generating a 
geometric model of the core stage and SRB. Since CGT reads files 
in PLOT3D format, a simple text file was created that contained the 
information for a single point located at the origin. This file was 
named origin.dat and details for this file are in the appendix A 
(figure Al). Once origin.dat was created, the file was opened in 
CGT’s Graphical User Interface (GUI) named OVERGRID. There 
was a read append feature in OVERGRID that allowed additional 
files to be appended to the current OVERGRID session. Thus, the 
file origin.dat was appended several times and each appended point 
was then translated from the origin to locations where there were 
geometric discontinuities in the rocket’s geometry. Since the core and booster were axisymmetric bodies, one 
feature curve was needed to create the core and SRB respectively. Figure 2 is an illustration of the surface grid 
generation process employed to create the surface grid for the core stage. A similar process was used to create the 
SRB surface grid. 

Due to NE-M1 unfamiliarity with the OVERFLOW software and its dynamic mode, a half-configuration of the 
SRB simulation was first created to experiment with parameters such as applied forces and moments, and time- 
steps. A full simulation that included both SRBs was constructed after the correct OVERFLOW files were created. 
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Table 1. Initial flight conditions used for 
the SRB separations simulation and 
LPVD simulations. 



Figure 2. The surface grid generation process of the core stage in CGT: a) Feature points, b) Feature curve, 
c) Feature curve revolved about x-axis. d) Points added and redistributed on surface. 


B. Geometric Modeling of Launch Pad Vehicle Drift Simulation 

Since the project was intended for demonstrational purposes, a decision was considered to make LPVD 
simulation two dimensional (2D). The launch pad that was modeled was derived from a mobile launch platform that 
was designed for two SRBs and three RS-25E rocket engines. Figure 3(a) illustrates an isometric view of the launch 
pad and figure 3(b) shows the 2D view of the launch pad as was desired for modeling. 
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For the CFD model of the launch pad, only the bulk dimensions, 
such as the total heights and widths were implemented. The 
intricacies of the truss structure were not represented in the CFD 
model. Moreover, the CFD model of the launch pad does not allow 
air to flow through the tower as it actually does. 

Since our proposes were to model the SLS under the least 
favorable wind conditions and due to the way the launch pad was 
designed to hold the launch vehicle, reducing the problem to two 
dimensions required the removal of the SLS’s SRBs from the CFD 
model. The least favorable wind condition would be a twenty knot 
crosswind that causes the launch vehicle to drift toward the tower. 

The geometry for the pad was created in OVERGRID in a similar 
manner as the geometry from the SRB separation simulation. Instead 
of using the origin. dat file that was created for the SRB separation 
simulation, a new PLOT3D formatted file was created and it was 
named 12pts.dat. This file contained the twelve points needed to 
represent every corner of the launch pad model. See the appendix A to 
view the contents of I2pts.dat. Figure 4 outlines the grid generation 
process for the launch pad. Note the x, y, and z coordinate directions 
in Figure 4(d). On the other hand, the geometry of the core was 
derived from the same grids that were developed for core stage in the 
SRB separation simulation. As shown in figure 5(b), the first and last 
grid planes were extracted from the volume grid of the core. Then 
those two planes were concatenated to form a single grid as shown in 
figure 5(c). Finally, the grid was translated and rotated to sit on the launch pad grids as shown in figure 5(d). Note 
that the origin is on the centerline of the core stage and at the top surface of the launch pad. 
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Figure 3. The launch pad that was 
modeled for the LPVD simulation. 

(a) Isometric view of launch pad. b) 2D 
view of tower for CFD model. 



Figure 4. The grid generation process for the launch pad CFD model as viewed in OVERGRID, a) Points in 
PLOT3D fde named 12pts.dat. b) Two feature curves created from concatenating points together, c) Points 
added between initial points for finer grid spacing, d) Resulting near-body grids for launch pad. 
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a) b) c) d) 

Figure 5. Grid generation process for the core stage of the 2D LPVD case, a) Near-body volume grid 
generated for the 3D SRB separation case, b) First and last grid planes extracted, c) Concatenated grid 
planes to form near-body volume grid, d) Near-body grid rotated, translated, and read into launch pad grid 
file. 


III. Requirements for Simulations in OVERFLOW Mode-D with Grid Movement 

OVERFLOW makes use of Mode-D to enable grid movement; also known as OVERFLOW-D. This section will 
summarize the OVERFLOW-D input requirements and the procedures needed for both the SRB separation and 
LPVD simulations. Typically a new directory (referred to as working directory in this report) was created that would 
contain all of the required input files for each simulation case. 


A. Near-Body Grids Represent Solid Geometry 

OVERFLOW-D only requires the near-body grids of a problem to simulate RBD. Thus, a PLOT3D file named 
grid, in that contains the near-body grids for a particular simulation was required in the working directory. 

As previously mentioned, both a full- and half-configuration model was setup for the SRB separation simulation. 
Thus the grid. in for the half-configuration contained half of the core stage and a single SRB. The grid, in for the full 
case contained the full volume grid for the core-stage and two 
SRB near-body grids. Due to the ability to translate (and rotate) 
grids to other locations in the XML files, both volume grids for 
each SRB was placed in the same coordinate location. See the 
XML file section below. 

The grid, in for the LPVD simulation contained the near- 
body grids for the two launch pad components and the core. 

Note that since this case was 2D, OVERFLOW required a total 
of three grid planes to be supplied. Thus the GRIDED utility 
(batch mode) in CGT was used to add a plane in front and 
behind the original plane that was located at y = 0. The two 
additional grid planes supplied had to be located at y = 1 and 
y = - 1 respectively. Figure 6 shows the graphical contents the 
grid. in file for the 2D LPVD and highlights the three planes 
needed. 



Figure 6. Illustration of the three planes 
needed for a 2D simulation in OVERFLOW. 


B. Force and Moment Coefficient Pre-Processor 

MIXSUR is a CGT utility that was needed to generate input files that OVERFLOW-D utilized to compute the 
force and moment coefficients at each computational iteration, grid. in had to be in same directory as the directory 
that ran MIXSUR so that an integration surface could be generated. A file named mixsur.inp was needed for each 
simulation. Reference 5 was frequently consulted during the generation of the each mixsur.inp file. The same 
mixsur.inp that was created for the 3-DOF LPVD case was reused for the 2-DOF case. For the SRB separation 
simulations, one mixsur.inp was needed for the half-configuration. For the full-configuration case, the same 
mixsur.inp for the half-configuration case was modified to include the second SRB. See the appendix for the 
mixsur.inp files created for each case. 
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One of the many files that were generated after running the 
MIXSUR utility is grid.i.tri. For each simulation, this file was 
opened in the OVERGRID GUI and visually inspected before 
running OVERFLOW. Figure 7 is the grid.i.tri file that was 
generated for the half-configuration SRB separation case. 

C. OVERFLOW Input File 

For every simulation, OVERFLOW read a text file that 
contained all of the CFD input variables. The default name was 
over. namelist or overflow, inp. Different input files were 
distinguished by naming the files of the form overflow*. inp 
where * denotes text that distinguished one input file from 
another. This naming convention was very useful when a few 
changes, such as angle-of-attack and freestream Mach number 
were made, and also when a shell script was created to 
automate running multiple cases. 

The OVERGRID GUI in CGT allowed the overflow. inp 
file to be generated by interactively selecting the OVERFLOW flow solver input options. After an OVERFLOW 
input file was generated by OVERGRID, it was often the case that several parameters needed to be experimented 
with. 



Figure 7. Force and moment integration 
surface for the half-configuration SRB 
separation simulation that was created by the 
CGT utility named MIXSUR. 


D. Configuration and Scenario XML Files in OVERFLOW Mode-D 

There are three ways in OVERFLOW-D to 
specify dynamic grid motion 1 . The method chosen 
for this project was to set dynamical parameters in 
Extensible Markup Language (XML) files that use 
the Geometry Manipulation Protocol 6 (GMP). For 
each computation in OVERFLOW-D, the XML files 
named Config.xml and Scenario.xml had to be 
present in each working directory. For details of the 
contents in the XML files, refer to chapter 5 of 
reference 1 and for the XML files made for each 
simulation in this project see appendix B, C, and D. 

Grouping grid numbers to a component and any 
translations or rotations were contained in the 
Config.xml file. The Scenario.xml had the ability to 
specify inertial properties, aerodynamic 6-DOF, 
forces and moments, motion constraints, and any 
additional translations or rotations. Table 2 
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±513,200 
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±2,566,000 
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Table 2. Inertial parameters 
separation simulation. 


used in the SRB 


summarizes the inertial parameters used for the SRB in SRB separation simulations. Note that the F y and M z values 
in table 2 were positive or negative according to which SRB the load was being applied to. See chapter 5 of 
reference 1 for a list of equations that non-dimensionalize the dynamics quantities by the reference velocity. As a 
final comment for the non-dimensional values in table 2, the reference velocity used for all SRB separation 

simulations was 4,112 ft/s. The reference velocity chosen was the 
relative ground velocity that was stated in reference 4 however, the 
relative wind velocity of 4,106 ft/s in reference 3 should have been 
used instead. The OVERGRID GUI was used to assist in the 
construction of each XML file. 

Table 3 lists the inertial properties and their non-dimensionalized 
values that were input into the Scenario.xml file for the LPVD 
simulation. F z was a constant value that simulated the thrust that the 
launch vehicle would produce and was calculated by assuming a 
thrust to weight ratio of 1 .46. The reference velocity was determined 
Table 3. Inertial parameters used in the from the Mach "umber that would represent a 20 knot wind velocity 
LPVD simulation. and s P ee d of sound in air at sea level. 
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E. Object X-rays in OVERFLOW D-Mode and their Function 

Hole-cutting is a standard procedure that is needed in Chimera grid systems. Hole-cutting is the name of the 
procedure for blanking-out grid points. As several grids are overlapped, some grids may have points that are defined 
within solid geometry; such points must be removed or blanked-out of the computational domain. Figure 8 is an 
illustration of hole-cutting before and after the hole-cutting procedure and after grid motion has occurred. 
OVERFLOW-D has a built-in hole-cutter. 

The method in which holes were cut in OVERFLOW-D was by supplying object X-rays 1 . The X-rays that were 
required for both simulations were created in OVERGRID and the filename had to be named xrays.in. 

A CGT utility named XRAYED was a very noteworthy tool in the process of creating a moving-body simulation 
in OVERFLOW-D. Each near-body grid that defined a hole-cutter needed an X-ray and each X-ray had a 
component ID. The component ID was an integer that had to match the n ,h component defined in the Config.xml file 1 . 
It was crucial that the component ID matched the order of the defined components in the Config.xml file. For 
example, if the component ID of the third defined component in a Config.xml file is not equal to 3, then, if and when 
the grid(s) of the third component moves under the influence of dynamic motion, then the X-ray will not move with 
it; thus, the hole-cutting will be ineffective. Therefore, XRAYED was needed during the problem-stepup phase of 
each simulation to make sure that the component IDs of each X-ray were properly numbered. 

For the half-configuration SRB separation case, two X-rays were required; one for the core and one for the 
booster. For the full-configuration of the SRB separation case, XRAYED was used to duplicate the X-ray that was 
already defined for the booster in the half-configuration case and a new X-ray was needed for the core. For the 
LPVD case, one X-ray was used for the near-body grid of the core, and, since both near-body grids for the launch 
pad would remain fixed, one X-ray was used for both launch pad grids. 



a) b) c) 

Figure 8. Hole-cutting required for the LPVD simulation, a) Before hole-cutting procedure, 
b) After hole-cutting procedure, c) Hole-cut after near-body grid for the core stage has moved. 
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IV. Results 

Fully viscous solutions were computed for both the SRB separation and the LPVD simulations. For each 
simulation case, a steady-state solution was first computed until the aerodynamic coefficients stabilized. Then the 
solution and grid file was copied into a new working directory that contained all of the input files required for the 
dynamic simulation and the restart feature in OVERFLOW was activated. 

A. Half-Configuration of SRB Separation 

A steady-state solution was computed for the half-configuration SRB separation case. First a subsonic case was 
completed to ensure all OVERFLOW input parameters were selected correctly. When a supersonic case was first 
assembled and computed, the solution diverged. Thus, a smoothing parameter in the overflow. inp namelist file 
named DELTA was changed from the default value of 1.0 to 7.0. Figure 9(a) shows the subsonic solution at the 
plane of symmetry for the steady-state case. Figure 9(b) shows the solution at the desired Mach number of 3.877. 



Figure 9. Steady-state solution of the SRB separation case. Solutions are displayed at the symmetry plane, 
a) M a0 = 0.8, a = 0.0° b) = 3.877, a = 0.0° 


Interpolation error was evident in the wake region of the core in figure 9(b). This error was attributed to the 
difference in cell size form the core’s near-body grid and the nearest off-body grid. The near-body grid of the core 
shows a Mach number of about 3.2 and then an abrupt increase to about 3.6. 

With a time-step of 0.001 second, the simulation was allowed to run for a total of 5000 iterations which 
simulated 5 seconds. Figure 10 shows the predicted grid motion of the right SRB. Figure 1 1 shows the coordinate 
values of the SRB CM in all three dimensions. The relative orientation of the coordinates axes to the vehicle is also 
shown in figure 1 1 . The applied ejection forces and moment (table 2) were applied to CM from 0.0 to 0.75 seconds 3 . 
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XYZ Coordinates of SRB CM; SRB Separation; Half Configuration 
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Figure 11. Plot of the SRB CM in the x, y, and z directions during SRB separation for the half- 
configuration. At = 0.001s, M 00 = 3.877, a = 3.0°, 7 = 36.996°, Reoo= 18,370.77, p = 0.0° 
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Figure 12. Convergence of steady-state lift and drag coefficients for the full- 
configuration of the SRB separation case. 
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B. Full-Configuration of SRB Separation 

Once all issues were resolved in the half-configuration case for the SRB separation, a full simulation was 
assembled. Figure 12 (page 9) shows the convergence of the aerodynamic coefficients for the steady-state case. Due 
to the larger size of the grid file, the 6-DOF solution calculation was proceeding at an exceedingly slow rate of 1 ,000 
iterations every 2 days on 4 processors. A 100 processor computation was then submitted to the NASA Columbia 
Server at Ames Research Center. The computational time was significantly reduced to about 20 minutes for 1 ,000 
iterations. 

C. LPVD Simulation in 3-DOF 

The first simulation created for the LPVD was a 3-DOF case. The simulation results predicted an unfavorable 
motion as shown in figure 13. The motion was attributed to the low principal moment of inertia about the y-axis that 
was selected (see table 3). This simulation could be considered to be a rocket that launches form a launch pad 
without control systems. 



Figure 13. Velocity contours for the 3-DOF simulation of the SLS as it departs from a launch pad. A steady 
20 knot crosswind flows from left to right. Time is increasing from left to right then top to bottom. The 
motion was not expected. 


D. LPVD Simulation in 2-DOF 

Once the unexpected motion was viewed in the 3-DOF case, the Scenario.xml file was modified to constrain the 
rotation so that only 2-DOF remained (see appendix D, figure D4). As shown in figure 14, the predicted motion 
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seemed reasonable. The core stage drifted about 1.2 feet in the direction of the tower following 8 seconds of flight 
time. The thrust was not activated until 2 seconds of simulation time. The time step for this simulation was 0.001 
second. Figure 1 5 shows the flow solution for the 2-DOF case. 


X Coordinate of Core CM; LPVD; 2-DOF 



Z Coordinate of Core CM; LPVD; 2-DOF 

650 i 1 1 1 1 1 1 1 r 



123456789 10 11 

Time | si 

Figure 14. X and Z coordinates of launch vehicle CM for the 2-DOF case. Thrust was 
activated after 2 seconds of simulation time. (Note that the highest point on the launch 
tower was 351 feet. (See appendix A; figure A2) 
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Figure 15. Velocity contours for the 2-DOF simulation of the SLS as it departs from a launch pad. A steady 
20 knot crosswind flows from left to right. Time is increasing from left to right then top to bottom. 


V. Concluding Remarks 

NE-M1 has not had the capability to simulate RBD with interactions from accurate aerodynamic load models. 
During the course of this project, two simulations were prepared in OVERFLOW-D to demonstrate the influence of 
aerodynamic loads on rigid bodies in motion. Since a SLS is being considered by NASA, the geometry for the 
simulations were based on conceptual SLS designs. 

A fully 3D model with 6-DOF was developed to simulate the dynamical motion of SRBs as they separate from 
the core stage. Although the accuracy of the simulations were not verified, a template that contains lessons learned 
from using 6-DOF motion in OVERFLOW-D can be created so that future users can avoid issues encountered in 
this project. 

LPVD of the conceptual SLS was also simulated in 2D by allowing 2-DOF and 3-DOF. Since the objective of 
the project was for demonstrational purposes, the decision to reduce the simulation to 2D allowed for rapid 
preparation and computation of the simulation. 

The inertial properties employed in these simulations were not accurate; however, NE-M1 can now begin to 
approach more detailed CFD analyses that involve RBD. In the future, Guidance, Navigation, and Control experts 
should be consulted for an accurate description of the ejection loads that should be applied to the SRBs. 
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Appendix 


A. File contents of origin.dat and 12pts.dat 


Figure Al. PLOT3D formatted file 
named origin.dat for use in OVERGRID 
to create a single point at the origin. 


Figure A2. PLOT3D formatted file made for 
use in OVERGRID to create 2D launch pad 
model, (each line of text that has multiple 
significant figures are the x, y, and z 
coordinates of the corners in feet.) 
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47 . 75O000000E+00 
47 . 750000000E+00 
89 . 750000000E+00 
89 . 75O000000E+O0 
102 . 0O000000E+G0 
102 . 00O000O0E+00 
-56 . 0O00000E+0O 


-24 . 000000OE+00 
-24 . GG0GOG0E+GG 
0 . 0OOO0OOE+GO 
0 . 000O000E+0O 
35 1 . 0000000E+00 
35 1 . O00G000E+00 
75 . OO00O0OE+00 
0 . GGO00G0E+00 
-24 . 0G0O000E+G0 
24 . GG00G00E+O0 


111 

0.0 0.0 0.0 


B. OVERFLOW-D Input Files for Half-Configuration of the SRB Separation Case 


(GLOBAL 

RESTRT ■ T., NSTEPS - 5000, NFOMO - 1. WAVE - 1000, 
NOT - 102, NQC - 0. 

TPHYS - 0, DTPHYS- 4.112, 

HULTIG • T . FHG « F . FHGCYC - 150,1(0, 

(END 

(OMIGLB 

I RUN - 0. 

I600F - 2. OYNHCS - .T., N ADAPT - 200. 

(END 

(DCFGLB 

DOUAL - l 0, NOR FAN - 0, 

(END 

(GBR1CK 

0BGR1DS - T.. 

OFAR - 3000, DS - 0.9. 

XNCEN - 160, YNCEN - 00. ZNCEN • 0.0. 

1 YMIN - 1. P YMIN « 0.0. MIN8UF - 4, 

(END 

(BRK1NP 

N BRICK - 0. 

(END 

(GROUPS 

(END 

(XR1NF0 1DXRAY- 1. IGXLIST- 2. XDELTA - 1.25, (END 

(XRINFO IDXRAY- 2. IGXLIST® l. XDELTA - 1.25. (END 

(XRINFO 10XRAY- 1. IGXLIST- -l, XDELTA - 17, (END 

(XRINFO 1DXRAY- 2. IGXLIST- 1. XDELTA - 6, (END 


(FL01NP 

FSMACH - 3.877, REFNACH - 3.877. ALPHA - 3.0, 
REY - 18370.77. TINF - 468 2408. 

(END 


(VARGAM (END 


(GRONAM 

NAME - ‘Min'. 
(END 


(MEIPRM 

IAHS - 6. 1LHS - 6, IOISS - 3. 

BININ - 10. 

(ENO 

(TIMACU 

ITIME - 0, CFLMIN - 5.0, CFLHAX - 0.0. 

(END 

(SHOACU 

ISPEC - 2, SMOO • 10, D1S2 - 2.0. DIS4 - 0 04. 
DELTA - 7, 

(END 


(VISINP 

V1SC - T , 
(END 

(BCINP 
IBTYP « 
IBOIR - 
JBCS - 
3 BCE - 
KBCS « 
KBCE - 
LBCS - 
LBCE - 
(END 


(BCINP 
IBTYP . 
IBOIR ' 
JBCS ■ 
JBCE • 
KBCS - 
KBCE • 
LBCS • 
LBCE • 
(END 


(TIMACU 

(ENO 


(SHOACU 

(ENO 


(VISINP 

vise • 


17. 

2, 

1. 

•1. 

1. 

1. 

1. 

-1. 


(GRDNAN 

NAME - ‘booster* , 
(ENO 


(METPRM 

(ENO 


NAME ° Off-body grids'. 


(GADNAM 
NAME 
(END 

(NITERS (ENO 

(METPRM 

(ENO 

(TIMACU 

(ENO 


(VISINP 

vise ■ 

(END 


Figure Bl. Overflow. inp for half-configuration of the SRB separation case. 
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|c?xml versions’ 1.0* encodings’ utf 8' ?> 
Configuration AngleUnits H degree"> 

Component Names ‘main’ 1 Types"st ruc“> 

<Data> Grid List=l </Data> 

</Component> 

Component Na«e= ’booster" Types" struck 
<Data> Grid Lists2 </Data> 

</Component> 

</Configuration> 

Figure B2. Config.xml file for half- 
configuration of the SRB separation case. 


<?x«l versions' l . o ’ encodings' utf -8'?> 

Scenario AngleUnits "degree" Gravitys-l.43441E-6,0.0, -1.20217E-6"> 

<Aero6dof Co«ponent="boostef " Start="0.O" Duration="0.0"> 

<InertialProperties 
Nasss‘3 . 7763e+09“ 

CenterOfHasss" 169. 484b, 23.4667, 0.0" 

PrincipalNoaents0fInertia="143 191840543.0, 7510971307750.0, 7522448851710.0“ > 

<PrincipalAxesOrientation Axis=" 0.999999977151, -9.45999818456e-5, 0.00019169779858* Angle="0. 0"/> 
</InertialProperties> 

<AppliedLoad 

Start="0.0" Durations"750" Fraae="body“ Forces“l 6e4, 2.0e4, 0.0" Hoaent="0.0, 0.0, -1.0E5" 
</AppliedLoad> 

</Aero6dof> 

</Scenario> 

Figure B3. Scenario.xml file for half-configuration of the SRB separation case. 


e, 

586, 506, 

• 1 , 

, 0, 6 

FSHACH , ALPHA, BETA, REY , GAMINF , TINF 

2 





NREF 

33 

, 855.3, 160, 

6 ., 0 . 

REFL,REFA,XHC,YMC,ZMC 

12 

.2, 

116.9, 

169.4845, 23.4667, 

6.0 REFL , REFA , XMC , YMC , ZHC 

2 





NSURF 

1 , 

1 




NSUB, IREFS 

1 , 

3, 

1,-1, 

1 , 

- 1 , 1 , 1 

NG,IBDIR,JS,JE,KS,KE,LS,LE 

0 





NPRI 

1 , 

2 




NSUB, IREFS 

2, 

3, 

1,-1, 

1 , 

- 1 , 1 , 1 

NG,1BDIR,JS,JE,KS,KE,LS,LE 

0 





NPRI 


2 

NCOMP 

nain 

1, 1 

NIS, IREFC 

1 

(surface nuBbers) 

booster 

1, 2 

NIS, IREFC 

2 

(surface nuBbers) 


Figure B4. Mixsur.inp file for half-configuration of 
the SRB separation case. 


C. OVERFLOW-D Input Files for Full-Configuration of the SRB Separation Case 
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Figure Cl. Overflow.inp file for full-configuration of the SRB separation case. 


<?xml versions' 1 O' encodings' utl 8' ?> 

<Configuration AngleUnits"degree“> 

Component Na«es"main" Types"struc“> 

<Data> Grid Listsl </Data> 

</Component> 

Component Na*es"Rbooster " Types’st ruc"> 

<Data> Grid List=2 </Data> 

</Component> 

Component Na«es"Lbooster'' Types-struc^ 

<Data> Grid List=3 </Data> 

<T ransform> 

translate Displacements"© .0, 46.9338, 0.G" /> 

</Transform> 

</Component> 

</Configuration> 

Figure C2. Config.xml file for full-configuration of the SRB 
separation case. 


|c?xml versions' 1.0' encodings' utf -8' ?> 

<Scenario AngleUnits "degree" Gravity="1.43441E-6,0.O, 1.20217E-6"> 

<Aero6dof Conpor»ent="Rbooster " Starts"©. 0" Durations"© . 0*> 
dnertialProperties 
Hass=“3 . 7763e+09" 

CenterOf Hass =*169.4845, 23.4667, 0.6“ 

PrincipalMonentsOf Inertias ”143 191840543 .0, 7510971307750.0, 7522448851710.0" > 

<PrincipalAxesOrientation Axis=’-0. 999999977151, -9.45999818456e-5, 0.00019169779858" Angle="0.G"/> 

</ Inert ialProperties> 

<AppliedLoad 

Starts "0.0" Durations'' 750" Fra»e="body " Force="1.6e4, 2.0e4, 0.0" Ho«ent="0.0, 0.0, -1.6E5" 
</AppliedLoad> 

</Aero6dof> 

<Aero6dot Coaponents'Lbooster" Starts’©. 0" Durations "0 . 0“> 

<InertialProperties 
Mass=”3 . 7 763e+09 " 

CenterOf Hass = " 169 . 4845 , -23.4667, 0.0" 

PrincipalHoaentsOf Inertias "143 191840543.0, 7510971307750.0, 7522448851710.0" > 

<PrincipalAxesOrientation Axis=" 0.999999977151, -9.45999818456e 5, 0.00019169779858" Angle="0.0"/> 

</ Inert lalProperties> 

<AppliedLoad 

Start=‘0.0" Duration="750" Fraae=’body'' Force="1.6e4, -2.0e4, 0.0“ Hoaent="0.0, 0.0, 1.0E5* 
</AppliedLoad> 

</Aero6dof> 

</Scenario> 

Figure C3. Scenario. xml file for full-configuration of the SRB separation case. 
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e. 506, 500, -1, 0, 0 
3 

33, 855.3, 160, ©., 0. 

12.2, 116.9, 169.4845, 23.4667, 
12.2, 116.9, 169.4845, 23.4667, 
3 

FSMACH , ALPHA , BETA , REY , GAMINE , T INF 
NREF 

REFL,REFA,XMC,YMC,ZMC 
0.0 REEL , REFA, XMC , YHC , ZMC 
0.6 REF L, REF A, XMC, YMC, ZMC 
NSURF 

1, 1 

1. 3, 1,1, 

e 

1, -1, 1, 1 

NSUB, IREFS 

NG, IBDIR, JS, JE ,KS ,KE, LS, IE 
NPRI 

1, 2 

2. 3, 1,-1, 

0 

1, -1, 1, 1 

NSUB, IREFS 

NG, IBDIR, JS,JE,KS,KE,LS, LE 
NPRI 

1. 3 

2, 3. 1,-1, 

0 

1. -1, 1, 1 

NSUB, IREFS 

NG. IBDIR. JS.JE.KS.KE.LS.LE 
NPRI 

3 


NCOMP 

Min 
1. 1 
1 


NIS.IREFC 
(surface numbers) 

R booster 
1, 2 
2 


NIS.IREFC 
(surface numbers) 

L booster 
1. 3 
3 


NIS.IREFC 
(surface numbers) 


Figure C4. Mixsur.inp file for full-configuration 
of the SRB separation case. 


D. OVERFLOW-D Input Files for the LPVD Case 



Figure 1)1. Overflow.inp file for both 2-DOF and 3-I)OF LPVD case. 
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<?xml versions’ 1.0’ encodings' utf -8' ?> 
<Configu ration AngleUnit="radian"> 

cComponent Name="core" Type="st ruc M > 
<Data> Grid List=l </Data> 
</Component> 

cComponent Name= ,, pad'’ Type=”st ruc"> 
<Data> Grid List=2, 3 </Data> 
</Component> 

</Configuration> 

Figure D2. Config.xml file for both 2-DOF and 3- 
DOF LPVD case. 


<?xml versions' 1 . 0 * encodings' ut 18 ‘ ?> 

<Scenario AngleUnit="radian" Gravity=”0 . 0, 0.0, -0 . 02823566 "> 

<Aero6dof Component="core" Start="0" Duration="0"> 
dnertialProperties 
Hass="7 .53E7" 

Center0fMass="0, 0, 160" 

PrincipalMomentsOf Inertias"© .0, 7E7, 0.0" > 

<PrincipalAxesOrientation Axis="1.0, 0.0, 0.0" Angle= ,, G.0"/> 
</InertialProperties> 

<AppliedLoad 

Starts"©.©" Frames" body" Forces"©.©, 0.0, 3.1E6" Moments”© .0, o.o, o.o" 
</ApplledLoad> 

</Aero6dof> 

</Scenario> 

Figure D3. Scenario.xml file for 3-DOF LPVD case. 


<?xml versions' 1.0* encodings' utf -8' ?> 

cScenarlo AngleUnits" radian" Gravitys"0, 0, -0 . 0282357 "> 

<Aero6dof Components" core" Starts"© " > 
dnertialProperties 
Masss"7 . 53E7 " 

CenterOfMass="0, O, 160" 

PrincipalMomentsOf Inertias"©. 0, 7E7, 0.0" > 
cPrinclpalAxesOrientatlon Axiss"l.G, 0.0, O.O" Angles"©. 0“/> 

</ Inertia IP rope rties> 

<ApplledLoad Start="G.O " Durations" o.o - Frames "body" 

Forces" 0.0 , 0.0 , 3.1E6 " Moments" 0.0 , 0.0 , 0.0* /> 

constraint Starts"©. 0" Du rations “0 . 0" Translates"©, 0, 0" Rotates"©, 1, 0" /> 
</Aero6dof> 

</Scenario> 

Figure D4. Scenario.xml file for 2-DOF LPVD case. 


0, 500, 500, -1, 0, 0 FSMACH, ALPHA, BETA, REY.GAMINF.T INF 

1 NREF 

33, 10378.5, 0.0, 0.0, 160. REFL.REFA, XMC , YMC , ZMC 

1 NSURF 


1, 1 NSUB ( core) , IREFS 

1, 3, 1, -1, 1, -1, 1, 1 NG, IBDIR, JS, JE,KS,KE, LS, LE 

0 NPRI 


1 NCOMP 


core 

1, 1 NIS.IREFC 

1 (surface numbers) 


Figure 1)5. IMixsur.inp file for both 2-DOF and 3-DOF 
LPVD case. 
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